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Based on the theory of Matrix Product States, we give precise statements and complete analytical 
proofs of the following claim: a large fractionalization in the magnetization or the need of long-range 
interactions imply large entanglement in the state of a quantum spin chain. 
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INTRODUCTION 

Entanglement plays a central role in many-body quan- 
tum systems as it can be used to understand the struc- 
ture of the quantum states that appear in nature. In 
systems governed by short-range interactions, low en- 
ergy states possess very little entanglement. In contrast, 
states evolved after quenches display large amounts of 
entanglement. These different behaviors, which are sup- 
ported by abundant numerical evidence, have been re- 
cently established on solid grounds in one spatial dimen- 
sion [1]. In particular, ground states of gapped (critical) 
Hamiltonians fulfill an area law, in which the entangle- 
ment entropy of any connected region is bounded by a 
constant (diverges at most like the logarithm of the num- 
ber of spins in that region) [2, 3]. These results immedi- 
ately imply that the ground state of a spin chain can be 
well approximated by Matrix Product States [4, 5], and, 
thus, such family of states captures the physics in one 
dimension [6-8]. 

Apart from the cases mentioned above, there exists 
practically no other physical situation where the exis- 
tence of large or small amounts of entanglement can be 
rigorously established. In this paper, we identify two sce- 
narios in one spatial dimension that can be connected to 
the presence of entanglement. We will restrict ourselves 
to systems described by MPS, and thus, our results do 
not apply to general situations. Nevertheless, given the 
fact that such family of states approximates well one di- 
mensional systems, we conjecture that our results are 
true in more general settings. 

The first scenario corresponds to the presence of frac- 
tionalization, a striking phenomenon that arises when- 
ever certain observables, which are expected to take in- 
teger expectation values, appear to be fractional-valued 
instead. The most prominent example of such behavior 
is the celebrated fractional quantum Hall effect [9-11]. 
In recent years, this phenomenon has been extensively 
studied in many systems, including spin chains [12, 13], 
where the magnetization per particle is fractionalized as 
a function of the external magnetic field. In the first 
part of this manuscript, we establish a lower bound for 
the entanglement entropy of any (connected and suffi- 
ciently large) region of a quantum spin chain in terms of 



the fractionalized magnetization. 

The second scenario corresponds to a situation where 
the area law does not apply, namely when studying a spin 
chain with long-range interactions. Intuitively, one can 
expect that such interactions give rise to large amounts of 
entanglement since any specific region will be correlated 
to any other part of the chain. However, it is very sub- 
tle to transform this intuitive idea into a rigorous result. 
The main reason is that the ground state of Hamiltonian 
containing long-range interactions may coincide with (or 
be very similar to) the ground state of another Hamilto- 
nian containing short-range ones and, therefore, fulfilling 
the area law. For instance, if we have an Ising model with 
decaying interactions (in the absence of a transverse mag- 
netic field), the ground state will be still a product state, 
which, in turn, is also the ground state of the Ising model 
with nearest-neighbor couplings. Such state does not dis- 
play any entanglement at all. Hence, we can only expect 
to have large amounts of entanglement whenever such 
examples do not exist; that is, whenever our state is (in 
some sense) not close to any other state corresponding to 
the ground state of a Hamiltonian with short-range inter- 
actions. In fact, we will prove a Theorem that formalizes 
this statement in the second part of this paper. 

In order to rigorously prove our statements, we will 
have to further develop the theory of MPS, extending 
previous results presented in Refs. [4, 14, 15], and deriv- 
ing new ones. Some of them are very intuitive, although 
the rigorous proofs are somewhat complicated. We will 
present in the main text of this paper the main steps and 
their intuitive interpretation, and leave for the appen- 
dices the technical details. 



MATRIX PRODUCT STATES 

This family of states describes a chain of N spins J, 
with d = 2 J + 1, and can be written as 
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E 

ii ,...,ijf=l 



tr [A n [!]••• A tN [N]} |*i - - - i N ) (1) 



Here, Ai[n] are D x D matrices associated to the spin 
in the n-th site of the chain. Our results, unless specifi- 
cally mentioned, concern translationally invariant states, 
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where Ai[n] = Ai independently of the site n. We will 
call the corresponding state \i^a)- 

Let us recall some known properties of such states 
(see, for instance, Refs. [4, 5]). MPS can be clas- 
sified into injective and non-injective. An MPS is 
called injective if there exists an L such that for re- 
gions of size L or larger, different boundary conditions 
give rise to different states; that is, the map T(X) = 

J2n i L tr [XA^ ■ ■ ■ A iL \ \ii ■ ■ ■ i L ) is injective. This is 

known to be equivalent to the fact that, after a suitable 
transformation of the form Ai XAiX^ 1 , one obtains a 
canonical form fulfilling the following conditions (that we 
will always assume for injective MPS): (i) ^ i AiA\ = 1, 

(ii) J^i A\h-AAi = Aa for a diagonal positive full rank 
matrix Aa, and (iii) the cp map defined as 

^ A {X)=Y J AiXA\ (2) 

i 

has 1 as its unique non-degenerate eigenvalue of maximal 
modulus. This canonical form is unique in the following 
sense: if A and B are matrices giving rise to different 
canonical representations of the same MPS, then they 
must be related by a unitary U according to e 10 Ai = 
UBiU^ . Non-injective MPS also possess a canonical form 
where the matrices are block-diagonal and the cp map 
associated to each block verifies conditions (i), (ii) and 

(iii) above, except for the existence of other eigenvalues 
of modulus 1. 



LARGE FRACTIONAL MAGNETIZATION 
IMPLIES LARGE ENTANGLEMENT 

Fractional magnetization in a spin chain occurs when- 
ever we have U(l) symmetry (generated, in this case, by 
an operator J z ), and the expectation value of the gener- 
ator m = {J z )/N, the magnetization per particle, fulfills 
that J — m = q/p, where p and q are coprime. When 
we change some external parameter, such as a magnetic 
field, the value of to generally changes in discrete steps, 
giving rise to typical plateaus in the magnetization. Our 
aim is to show that whenever a translationally and U(l) 
invariant MPS displays this phenomenon, the entangle- 
ment entropy of any sufficiently large region is greater 
than log(p). That is, the value of p imposes some lower 
bound on the entanglement present in the system. We 
will start out with a trivial example that will help us 
build an intuition about this statement, and then we will 
generalize this claim to arbitrary states. 

Let us consider J — 1/2, any two numbers q, p coprime, 
and construct a state of N = np (n integer) spins as 
follows. We consider first a p-particle state of the form 
| a) = | tt " " " III ' ' ' I): where q is the number of spins 
down. Then, we take n = N/p copies of such state, and 
build an equal superposition of the p possible different 



translations of \a}® n , 

w = -4 E r "»^ ( 3 ) 

VP m=0 

where r is the translation operator. This state is trans- 
lationally invariant, has U(l) symmetry generated by 
J z = J2 s n> where s z is the single-spin operator 2s z = 
I t)(t I -\l)(l I, i-e. JM) = N(l/2 - q/p)\ip), and, 
thus, exhibits fractional magnetization. Following the 
prescriptions of Oshikawa et al. [12], this example con- 
tains "periodic components" in order to display such a 
phenomenon. As one can see by simple inspection, if we 
take any region A of size L = kp, with k G {1, . . . ,n}, 
the reduced density operator can be written as 

1 

PL = ~ X) kmXVml (4) 
^ m=0 

where (p m are mutually orthogonal. Thus, the entropy 
of pl (and, consequently, the entanglement entropy) be- 
tween the region A and the rest is log(p). This toy model 
presenting such entropy is connected to the fact that, 
in this case, fractional magnetization arises because the 
ground state is a linear superposition of p-particle states 
which are both locally orthogonal (i.e. fully distinguish- 
able) and related through a translation. 

In what follows, we will consider the richer family of 
MPS in order to prove a related result. Note that the 
previous toy example is contained in the family of MPS 
just by considering the matrices 

^ = x>><i+ii, a t = E (5) 

i— 1 i—q+1 

However, general cases of MPS possess several difficul- 
ties. In particular: (i) finding a characterization of all 
MPS displaying fractional magnetization; (ii) the fact 
that, even if an MPS is a superposition of states related 
by a translation, nothing ensures that the reduced states 
will be sums of few pure and mutually orthogonal states. 
In any case, we are able to prove the following Theorem: 

Theorem 1. Let be a translational and U(l) invari- 
ant MPS of spin J, with magnetization per particle m and 
verifying J — m = | (p and q coprime). Then there ex- 
ists a constant 7 € N such that the entropy of the reduced 
density matrix of any region of size L = kjp (Vfc € N) 
verifies S(pl) > log(p), up to a exponentially small cor- 
rection in N — L and in k. 

In order to prove this, we proceed as follows: 

Proof. Let \ip) be an MPS, which is translational and 
U(l) invariant. We also impose that this MPS has spin 
J and magnetization per particle to, verifying J — to = ^ 
(p and q coprime) and consider its canonical form. If it 
has a single block, due to Lemma 8 in Appendix A, it 
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must be 7p-periodic, where 7 G N. This means that all 
the eigenvalues of magnitude one corresponding to the 
cp-map are the jp- roots of unit. Consequently, if we 
block -fp spins, then we can write the new matrices Ai 
as block-diagonal, with each block being injective and 
different (see Lemma 5 in Appendix A). We have now 
that the state can be written as linear combination 
(with equal coefficients) of jp different injective MPS, 
each of them being a translation of each other. In Lemma 
3 of Appendix A, we show that different injective MPS 
are orthogonal in the thermodynamic limit. Let L = k-yp 
(k G N), using Jensen's inequality we have that S(pl) > 
— log(tr(p| / )) which, by Lemma 4 in Appendix A implies 
S(pl) > log(7p) > log(p), up to an exponentially small 
correction in N — L and k, as in the example proposed 
above. 

If the MPS has many blocks in its canonical form, we 
will show that one can treat each of these blocks as in the 
single block case, obtaining an extension of the last result. 
Lemma 8 gives us 7 <E N such that all the blocks of the 
canonical form of \ip) have period jp. Let L = kjp, where 
k G N. We observe that the reduced density operator of 
a region comprising L sites verifies 

PL = ®i =1 l*pi, (6) 

up to a correction exponentially small in N — L and k 
(see Lemma 4 in Appendix A). The p,'s are the reduced 
density matrices corresponding to single blocks, where 
repeated blocks are simply reflected in the pi's. Using the 
single block case, we can ensure again that S(pi) > log(p) 
for all i. It is clear, from (6) and the subaditivity of the 
entropy, that S(pl) > log(p) up to another exponentially 
small correction in N — L and k, yielding the desired 
result. □ 

To prove the crucial Lemma 8, it will be enough to 
consider the characterization of symmetries for injective 
MPS [14], as well as an extension of the Lieb-Schutz- 
Mattis theorem for U(l) symmetry which is explained 
in Lemma 6. The first [14] will allow us to assert that 
all injective MPS corresponding to blocks must have the 
same symmetry and the same magnetization m. The 
later, that all these blocks should have a period multiple 
of p. Moreover, Lemma 5 in Appendix A ensures that 
states corresponding to different blocks are necessarily 
different. 



LARGE INTERACTION LENGTH IMPLIES 
LARGE ENTANGLEMENT 

Now, we turn to the other situation where one can 
prove the appearance of entanglement. For that, let 
us consider again a translationally invariant MPS, \iPa), 
which is not the ground state of any short-range (gapped 
and frustration-free) Hamiltonian. Furthermore, let us 
assume that it is also far away (as specified below) from 



any other state with this property for any given inter- 
action length. We will show that, as a consequence, its 
entanglement entropy will be large and, indeed, will scale 
with the range of the interaction. 

In fact, if we denote by p L A the reduced density operator 
of \ip) for a (connected) region containing L spins, we can 
prove the following Theorem 

Theorem 2. Let \4>a) be an MPS such that every 
state which is the unique ground state of a gapped 
frustration-free Hamiltonian with interaction length L 
verifies \\p^ — p L \\i > e. Then, for sufficiently large re- 
gions R, we have that the a-Renyi entropy S a (p A ) > 
aL + b log e + c, for a < g and where a, 6, c are constants 
depending only on the local physical dimension d of \ipA)- 

This claim can be proved by contradiction. We will 
suppose that for every connected region and for a < g 
we have that the a-Renyi entropy is upper bounded by 
an expression of the form aL + b\oge + c, for a < i, where 
a, 6, c are constants depending on the physical dimension 
of \iI>a}- It will be enough to prove that this implies the 
existence of a state, the unique ground state of a gapped 
frustration-free Hamiltonian with interaction length L, 
such that \\pa — p L \\i < e. 

The hypothesis on S a being small implies that we can 
find another MPS with a sufficiently small bond dimen- 
sion, D (in particular, D < d^^ 1 ^ 2 ) that is close enough 
to the original one. In order to do this, we will rely on 
[16, Lemma 2] and on a new bound for reducing the bond 
dimension of an MPS. More precisely, this bound will be 
of the form 

\\Pa - Pjh < 2\/2d i / 2 \/Z ( 5 1 / 4 + (2L + 3)6, 

where p L ^ will be the reduced density matrix which can be 

constructed from p\ by substituting the Kraus operators 

Ai by PA t P and A = PAP, where P = Y? It will 

be explained in further detail in Appendix B. 

Now, arbitrarily close to the MPS associated to the 
reduced density matrix there exists another which is 
the unique ground state of a Hamiltonian with gap and 
interaction length L. Taking into account that the inter- 
action length is closely related to the bond dimension at 
which the MPS reaches injectivity, this will be deduced 
from proving that all MPS (except for a set with measure 
zero) reach injectivity fast enough. Standard Algebraic 
Geometry, as explained in Lemma 11 of Appendix C and 
[17, 18], reduces this problem to finding the existence 
of a single MPS displaying this property. The existence 
of such an example can be obtained using quantum ex- 
panders, as explained also in Appendix C. 

A more detailed proof can be given as follows: 

Proof. Let us call Ai the ordered eigenvalues of p A , which 
can be taken as close as wanted to those of A ® A by 
enlarging region R [16, Lemma 2]. In this case, it is not 
difficult to see that, if we call Pi the ordered elements of 
A, then J2Zd+i P-i < J2Zd+i A * = : 6 - 
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Suppose that, for a = i and for all R, we can upper- 
bound the a-Rcnyi entropy by 

S a (p2) < | log e + 1 (L log d - log L) - log | (7) 

In Appendix C, we show that we can always construct 
a state, that we shall call \tp), of the form 



H,...,IL 
lL+l,---,iN 



( 8 ) 

where A i7 Bj,C k & M^ x ^, Ai = PAiP (being Ai the 
Kraus operators defining the original MPS), with bond 
dimension D < S L ~^I 2 and such that the fixed point 
for the associated channel is A = PAP, where we are 

considering that P — J^f 

In Appendix C, we also prove that all states of this 
form (except a set of measure zero) reach injectivity in 
L — 1 sites. Therefore, the one we have constructed in (8) 
is the unique ground state of a frustration-free Hamilto- 
nian with interaction length L [4, 5]. Using a straight- 
forward adaptation of [5, Section 6], this Hamiltonian is 
also gapped. Even though this state is not exactly trans- 
lational invariant, it verifies that its normalized reduced 
density matrix for particles 1 ... L is of the form 



J2 (^(a|[i il ...i i ,Ait i 

ii, ...,*!, \a,/3 
ji, — Jl 



...A] 



up to a exponentially small correction (see Appendix D). 

This will allow us to use a bound, which is proved in 
Appendix B, which states that 

\\Pa ~Pa\U< 2V2d L / 2 VLS 1 / i + (2L + 3)6 

< AV^d^VIS 1 /' =: e', 

since the first term in the sum is clearly larger than the 
second. It only remains to show that e' < e, or equiva- 

4 

lently, that 5 < 2 la l' 1L -jL ' ^ mce we nave taken R large 
enough, then up to a exponentially small correction in R, 
we can state that 

log(<5) < — ( SM) - log ° X 



a 



1 - a 



Using this and the fact that that D > S L ^I' 1 , it is 
enough to prove 

„ , „. 4a , logd /, 4a \ 

SM< T ^**e+-l-(l-—)L 

= ^loge+ -^-(Llogd- logL) - log | - log ^ 
5 11) 4 



where, in the last step, we have set a = g. This is given 
by hypothesis in Eq. (7) [19]. Therefore, there exists a 
state 1^), which is the unique ground state of a gapped 
frustration-free Hamiltonian with interaction length L, 
such that \\p£ — p L \\i < e, as we wanted to prove. □ 



CONCLUSION 

In this work, we have shown how MPS are power- 
ful enough to provide formal proofs of certain believed 
statements on strongly correlated spin systems that were 
lacking a mathematical treatment. In particular, we have 
stated and proved that, for the state of a quantum spin 
chain, either a large fractionalization in the magnetiza- 
tion or the impossibility of being well approximated by 
the ground state of a local Hamiltonian demands large 
entanglement. Moreover, since MPS seem to be the right 
representation for the low energy sector of ID systems, 
one may postulate the results being true in full generality. 
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Appendix A: Technical lemmas for the proof of 
Theorem 1 



Our first aim is to state and prove a couple of lem- 
mas formalizing the claim: "for injective MPS, different 
means orthogonal" . 

Lemma 3. Given two injective MPS, \ipA) and \iPb), 
then IHV'b)! = 1 up to an exponentially (in N) 

small correction. Moreover, either both are equal for all 
N, or trN-i J \'ipA)('>pB\ — up to an exponentially (in 
N — L) small correction. In particular, \(iPa\'4 , b)\ = up 
to an exponentially (in N) small correction. 

Proof. It is easy to see that {iPb\^a) = tr , where 

£a,b — J2i^i ® Si. Moreover, it is clear that the 
eigenvalues of £a,b are the same as those of the map 
Ka.b(X) = ^AiXBj, which gives = 1 

up to a exponentially small correction. To finish, it is 
enough to see that all eigenvalues A of Ka.b verify that 
|A| < 1. We will use the conditions verified by the canon- 
ical form of an injective MPS, that is: (i) X^i^^l = ^ 
(ii) J^i A\AAAi — Aa for a diagonal positive full rank 
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matrix A a, and (iii) the cp map Ea denned as 

E A (X) = J2A i XA\ (Al) 

i 

has 1 as its unique non-degenerate eigenvalue of maximal 
modulus. 

Let us take A such that ^2iAiXBj = XX, using (i) 
for A and (ii) for B we get 



|tr [XA B Xt] 



tr [AiXB\k B X^ 



< 



i V2 r 



E 1 



tr \xBjA B B l X' 1 ^ tr ^AA^A 1 A 



1 1/2 



= |tr[AA B At] |, (A2) 

where we have used the Cauchy-Schwarz inequality and 
tr [AAsA 1 '] > 0. So, if |A| > 1, we must have an equality 

and, therefore, aAflU, = A 1 J 2 B l X J< . Multiplying 
by the adjoint expression, summing in i, taking traces 
and using (i) and (ii) again we get that |a| = 1 and, 
hence, a = e %e . Finally, since As is invertible, we get 
Y,i BiX^XBj = At A, which, by (iii), leads to X*X = 1 
and implies that e ld Ai = XBiX^ . This means that \ip A ) 
and IV's) are equal up to a global phase, for all N. □ 



A similar proof gives the following 



Lemma 4. Given an MPS of the form |t/ 
S™=i 'MVv) such that the \ip r ) are different injective 
MPS, thentr [tf tf;] cx <5 r . s +0(e- L ) + 0( e -( JV - L )), being 
the reduced density matrix for L particles associated 
to |Vv)- 

The next thing we need is the following modification 
of [4, Theorem 5]. 

Lemma 5. Consider any MPS \4>a) & <C d ® N which has 
only one block in its canonical form with D x D matrices 
{Ai} and such thatEA has ft eigenvalues of modulus one. 
If (3 is a factor of N, then the state can be written as a 
superposition of (3 (3 -periodic different and injective MPS 
with equal coefficients and bonds Di (also with the prop- 
erty that Di = D). Otherwise, if (3 is not a factor of 
N, then \ip A ) = 0. 

Proof. The only thing to prove is that the (3 /3-periodic 
states are injective and different. In the proof of [4, The- 
orem 5] , based on [5] , one proves the existence of a set of 
orthogonal projectors {Pk} with ^ fc Pk — 1 such that 



E A (X) = '£p j E A (P j XP k )P k , 



(A3) 



E A has 1 with degeneracy exactly (3 as the only eigenvalue 
of modulus 1 , and each block in the block-diagonal form 
of the Kraus operators of E A given by (A3) corresponds 
to one of the /3-periodic states. Moreover, the space of 



fixed points is generated by Pk and the space of fixed 
points of the adjoint map is generated by P k AP k . 

The cp maps associated to the /3-periodic states are 
then E fe (A) = P k E^{P k XP k )P k (restricted to inputs 
with X = P k XP k ). It is clear that P k is its only fixed 
point, P k AP k the only fixed point of the adjoint map, and 
there is no other eigenvalue of modulus 1, which shows 
that all /3-periodic states are injective. Now, if two of 
them were equal, we would reach a contradiction in the 
following way. For simplicity, we reason in the case of 2 2- 
periodic states but the argument can be adapted straight- 
forwardly to the general case. E 2 A has block-diagonal 



Kraus operators of the form B t <g> 1 0) (0 1 + Cj (8) 1 1) (1 1 . By 
the hypotheses and the uniqueness of the canonical form 
for injective MPS, Bi = e ie UdW for all i. Then, apart 
from 1 <gi |0)(0| and 1 <g> we also get U <g> |0)(1| as 

an eigenvector of with eigenvalue of modulus 1; the 
desired contradiction. □ 

Finally, we need the following version of the Lieb- 
Schultz-Mattis theorem for U(l) symmetry. It is inter- 
esting to note that it does not use any MPS structure, 
so it is valid in full generality and for any spatial dimen- 
sion. Let us recall that, in [20, Lemma 17], we showed 
that any quantum state with a f/(l) symmetry given by 

the canonical generator of spin si' 1 ^ verifies that 



IV>> 



JgNm 



IV>> 



(A4) 



= with u g = e l9S *' ) and a magnetization per particle m. 



Lemma 6. Let m be any rational number and p E N 
such that there exist two quantum states of (local spin J 
and) pN and (N+l)p particles respectively, for some N, 
having both of them magnetization per particle m. Then 
p( J — m) = q with q integer. 

Proof. By expanding equation (A4) in the canonical 
basis, we get £ fel ... fc c kl ... kpN e l9 ^i k] |fci • • • k pN ) = 



J2 kl - 



^igpNm 



c kl ...k pN \ki ■ ■ ■ k pN ). Since it is a basis 



and the state is not zero, there must exist fci, ■ ■ ■ , k P M G 
{— J, — J + 1, ...J — 1,J} such that ^Zjkj = Npm. 
For the same reason, there must exist k[ , ■ ■ ■ , k' pN+p e 
{-J, -J+ 1, ... J- 1, J} such that J2j k 'j = (Np + p)m. 
Subtracting, we get that mp — J2j k'j — J2j % has the 
same character (integer or semi-integer) as pJ. □ 

With this at hand, if we consider an MPS \ip) of spin 
J and pN particles with a U(l) symmetry, given by the 

canonical generator of spin si' J \ we have the following 
lemma. 

Lemma 7. Let p be the smallest integer such that, after 
blocking p sites together, \ip) has a block- diagonal repre- 
sentation with injective blocks. Then p(J — m) = q, with 
q an integer. 

To see it we consider blocks of p-sites. From [14, The- 
orem 5], we know that each block is an injective MPS 
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with the same symmetry. Since, by Lemma 3, states 
corresponding to different blocks are equal or linearly in- 
dependent, all of them must have also magnetization m. 
Now, by the characterization of symmetries for injcctive 
MPS [14], we know that the matrices defining each block 
inherit the symmetry and therefore the associated MPS 
has magnetization m for all system sizes that are multiple 
of p. Lemma 6 finishes the argument. 
We also get a reciprocal. 

Lemma 8. Let us assume that J — m = | with 
gcd(p, q) = 1 in a U(l) symmetric MPS, then there ex- 
ists 7 € N such that the MPS has only jp-periodic blocks. 
Moreover (trivially from Lemma 5), states belonging to 
blocks of different periods are different. 

Proof. As above, all injective MPS corresponding to the 
blocks must have the same symmetry and the same mag- 
netization m. Therefore, Lemma 6 shows that only 
blocks of period multiple of p can appear. □ 

We are finally ready to prove Theorem 1: 

Proof. Let \tp) be an MPS, which is translational and 
U(l) invariant. We also impose that this MPS has spin J 
and magnetization per particle m, verifying J — m = ^ (p 
and q coprime) and consider its canonical form. Lemma 
8 gives us 7 <G N such that all the blocks of the canonical 
form of \tp) have period jp. Consequently, if we block jp 
spins, then we can write the new matrices A, as block- 
diagonal, with each block being injective and different 
(see Lemma 5). Using Lemma 3, the injective and differ- 
ent MPS associated to the blocks are also orthogonal in 
the thermodynamic limit. 

Let L = kjp, where k <G N. We observe that the re- 
duced density matrix of size L. verifies 



PL 



(A5) 



i=i 



up to a correction exponentially small in N — L, where 
we have used Lemma 3. Here, the pt's are the reduced 
density matrices corresponding to single blocks (giving 
rise to different states) and repeated blocks in this sum 
are simply reflected in the /Vs. 

Analizing the single block case, we can ensure that 
S(pi) > log(7p) > log(p) for all i, up to an exponentially 
small correction in N — L and in k. This is deduced by 
using Jensen's inequality so that S(pi) > — log(tr(/^)), 
and recalling Lemma 4. Using (A5) and the concavity 
of the Von Neumann entropy, it is clear that, if we have 
several blocks then S(pl) > min^ S(pi) > log(p) up to 
another exponentially small correction in N — L and k, 
yielding the desired result. □ 



Appendix B: Bounds on MPS approximation 

Let Ai e M.d be the canonical Kraus operators defin- 
ing an injective MPS, with A as its fixed point. We define 
the normalized reduced density matrix for L particles p\, 
up to a correction exponentially small in N — L, by 



p l a= E 



tr 



A] L ---A]KA t 



•■■AiL Ni ' ' ' i L){ji ■ ■ ■ 3l\ 



ii , 



(Bl) 

We will also define p 1 ^ as the normalized density matrix 
resulted of projecting the Kraus operators (and the fixed 
point) into a subspace of dimension D < D, that is, 

Ai = PA t P and A = PAP with P = Y%=\ l*><*l- E wil1 
be the cp map associated to A, and E the one associated 
to A. Taking all this into account, we can state and prove 
the following Theorem: 

Theorem. 

\\Pa - p\h < 2tr \K 1/2 ] VL6^ 4 + (2L + 3)5, 



\Pa~Pa\U < 2V2DVLS^ 4 + (2L + 3)5 



where 5 = tr 



A- A 



In order to do this, we must prove the following two 
Lemmas as preliminary results: 

Lemma 9. ||E L (A) - A||i < 2LS. In particular, 
tr [e l (A)1 > 1 - 2L5. 

Proof. Using both the definition of 6 and that E is con- 
tracture for the 1-norm, we get that ||A-E(PAP)||i < 5. 
The map P • P is also contractible for the 1-norm, so 

||A-PE(PAP)P||i 

< ||A-PAP||i + ||PAP-PE(PAP)P||i 

< 26 

This means that ||A - E(A)||i < 26, since E(A) = 
PE(PAP)P. However, E is also contractible respect to 
the 1-norm, so 

||A-E 2 (A)||i < ||A - E(A)||! + ||E(A) - E 2 (A)||i 
< 46 



The result can be obtained by induction. 



□ 



We will now define, under the previous notation for 
the Kraus operators and the fixed point, the following 
operators 



7 



a A = ^ tr 



4 



■ a\ aa u 

31 11 



\h ■ ■ -«l)0'i • • • Jl| 



n u 



<7 A ,p= Yl tr ^•••^ t 1 AA 4l ...A li P 



iii-jL 



|*i • ■ -«l)0'i ' ' ' jzj 



I 

where it is important to note that (74, p is a positive op- We obtain the result by collecting all bounds above. □ 
erator. 



Lemma 10. \\p A - p A h < \Wa,p - 4> A \\2 + + 
where (f>^ = tr E L (A) p A is the not normalized reduced 

density matrix generated by Ai . The same holds changing 
the 2-norm by the 1-norm in both sides of the inequality. 

Proof. By using the triangular inequality and the fact 
that || • ||2 < || • ||i, 

\\pA - P A \\2 < \\pA — &a\\i + \Wa - &A,p\\l 

+ \\<ta,p - <j>Ah + Ha - pa\U 

The first term can be calculated exactly 
\\pA -va\\i = tr A| L ---A| i (A-A)A il ---Ai ] 

il,...,ZL 

= s . 



Now 



\0-A,P - <Pa\\2 



< 



tr 
tr 



Q(£*) L {A® A)£ L Q 



tr 
tr 



Q{F*) L {A® A)F L Q 
Q(£*) L (A ® A)£ L Q 



Q{r*) L {k®K)T L Q 
where £ = Y,i A i®Au Q = P®Pand7*= (1®P)£(1 
P)- 

We have now all the necessary tools to prove the main 
Theorem: 



Proof of the Theorem. We start by bound- 



ing the term p 



tr 



Q(£*) (A <g> A)£ L Q 



This can be done by 



The first equality holds because the operator is positive 

and the 1-norm can be replaced by a trace and the second tr I Q(F*) (A ® A) T Q 

one holds because E is trace preserving. The second term adding and subtracting terms such that they differ in 
can be bounded in a similar way. 



one projector, i.e. 



Wa — ca.pIIi = tr 



P x Y, A\ L ...A\AA i ,...A iL P A 



L-l 



< S + tT 



p± E A\ L ...A\AA i2 ...A i 



V - S | tr [^ L Q(^"*) r_1 ^*(l ® P^){£*) L - r {A ® A) 

r=l 
L-l 

+ ^ Itr [f s (l ® P ± )£^ L - s - 1 0(£*) L (A ® A) 

s=l 



This holds because ||A - E(A)||i = ||E(A - A)||i = (5, 
since E is trace preserving and E(A) = A. Therefore, 
||o-a-cta,p||i <<5 + tr [P ± A] =26. 

Finally, the last term can be bounded using Lemma 9 Let us bound the first family of tcrms . B y 

applying the Schwarz inequality |tr AiBi] \ < 



because 



Ha-paIU 



-1 + tr 



E L (A) 



< 2SL 



tr 



itr = 



tr 



£ (VkA kl ...A kL PAl ■ • • AtAl ■ ■ ■ A] L _K^ ® A 1/4 ^ 



k 1 ,...,k L 
i\,...,i r 
yi, — ,3L- 



A 1 /*® AV4i fcl ■■■A kL A]. •••At .At P-LAt. ■••At. \/X 



8 



tr 



£ A 1/4 A jL _ r ■ ■ ■ A n A lr ■ • • A n PA[ L ■ ■ ■ A\kA kl ■ ■ ■ A kL PA\ x ■ ■ ■ A\a\ x ■ ■ ■ A] L _A^ ® A 1 



/2 



ki,...,k L 
il,...,i r 



tr 



Y, A 1 / 2 ® A 1/4 A fel • • • A fct 4 • • • ^_ i ^P x A] i • • • A] L _ r AA jL _ r ■ ■ ■ A J1 P ± A lr A lr _ 1 ■ ■ ■ A it A\ L ■ ■ ■ A\. A 



t A 1/4 



k lt ...,k L 
i\ ,. ..,i r 
Jl, — ,3L- 



The first term is equal to 

1/2 



tr 

< tr 



A 1/2 
Ai/2 



tr 



PE ij (A)PE L (A 1 / 2 ) 



1/2 



The second term is equal to 

1/2 



tr 



AV2 



tr 



E' r 



- 1 oE( 



_Ltt?L— r / 



P^E 



(A)P ± )t i (A 1 / 2 ) 



1/2 



< S^tv 



AV2 



where we have used that A < A (hence, 

tr P ± E i - r (A)P ± < S), and that both E and 

E are contractiblc for the trace norm. Therefore, 

2 



jj, r < tr 



A 1 / 2 



The result for the v s is exactly the 



Ai/2 



V6. 



same, so it follows that /i < 2Ltr 

The other term can be calculated in the same way, by 
replacing £ — > T and T — > £, and it gives exactly the 
same estimate. 

The second inequality follows from the first one, 



tr 



A 1 / 2 



< v Ptr [A], and the fact that ga,p — 4>a has 



□ 



rank < 2D, which then gives 

\\<ta,p - <t>Ah - ^2D~\\&A,P - ^lb- 
Appendix C: Injectivity can be reached fast 
We will prove here the following technical lemma 



Lemma 11. Every MPS (with the exception of a zero- 
measure set) of the form 

W= H tr(A il ---A iL B iL+1 C iL+2 ---C iN )\i 1 ---i N ) 

(CI) 

w/iere Ai,Bj,Ck € MdxD and L > 2 \°^ reaches injec- 
tivity in every region of length L — 1 . 



Proof. Since the set of MPS failing this property is clearly 
a projective algebraic subvariety of (fC D ® C D ® C d )® 3 , 
standard algebraic geometry tells us that, if this set is 
non-empty, since (C 13 <g> C D (8) C d ) 83 is irreducible then 
both projective varieties must be equal [21]. Therefore, 
it is enough to find a single MPS reaching injectivity as 
stated in this lemma, which has been verified numerically 
up to D = 200 and d = 50, and also analytically in the 
next Lemma of this Appendix using quantum expanders 
[22]. □ 

It is proven in [23] that for all d > 4, there exists a 
Hermitian trace-preserving completely positive map 



E(X) = A \ XA i 
i=i 

such that |A 2 | < (^^) (l + O (log(P)P^)), where 

Ai e M D . 

Take the MPS \ip) generated by the matrices Ai and 
consider the map 

V n (X)= tr i XA h ■ ■ ■ A in ] |*i • • • *„) 



We want to show 

Theorem 12. Assuming D is large enough, T n is an 
infective map for [24] 



n > 



fclog(P) 



+ 1, K = 8,d > 16 



log(d) 

This will be a consequence of the following 
Lemma 13. 

1 



sup 

tr[XtX]=l 



r n (x)tr n (x) - -tr [xtx] 



<P|A 2 |" 



Proof. Considering in Md the usual Hilbert- Schmidt 
Hilbert structure, it is easy to see that the LHS is equal 
to 



1 



r*r n 

D 



op 



mi 




Figure 1. Given the tensor A — ((a\Ai\/3))i a p which defines 
the MPS, and with the usual convention that rotating means 
complex conjugation, we can represent the map r*F„ as the 
map in the figure from systems ac to systems bd and the map 
E n as the same figure but now mapping systems cd to systems 
ab 



for the usual operator norm on the Hilbert space M. u ■ 

Moreover, in coordinates, calling £ = J2 i A i ^)A i , we 
have that 

r;r n - ll = ]T ( (cd\£ n \ab) - ^5 ab 5 cd } \bd)(ac\ . 

abed ^ ' 

just identifying M. D = C D ® C D and calling \ij) to the 
canonical (matrix) basis there. 

Since for each operator on an n dimensional Hilbert 
space, || • \\ p < || • 1 1 2 < \/n\\ ■ \\ op , being || • || 2 the Hilbert- 
Schmidt norm, and using that the Hilbert-Schmidt norm 
is invariant under arbitrary rearrangements of the coor- 
dinates, we get that 

1 



r*r 1 

D 



op 



< D 



1 



^ v cd\£ n \ab) - -S ab 5 cd )\ab)(cd\ 

bed 



l 



= £>||£"--|l)<l||| op 
= D||E" -E° 



= D\\£ n -£ c 
op = D\\ 2 \ n 



Op 



where we have used in the last step that E is hermitian 
and |1) denotes the unnormalized vector X)£=i 

□ 



Proof of Theorem 12. T n must be injective as long as 

1 



|A 2 |" < 



D 2 



(C2) 



Otherwise, taking a (normalized) X such that 
T n (X) = 0, we would get a contradiction to 
Lemma 13. Since we know from [23] that U2I < 



C 2 ^ 1 ) (l + O (log(D)D^^ it suffices to take n such 



that 



(Vp) (l + o(log(D) ^))" <5 L 



Taking logarithms 
21og(Z>)+nlog 
which is equivalent to 
n > - 



l + O log(.D)£>T5 



< 



21og(£>) 



log 



It is clear that taking D large enough we can upper- 
bound the RHS by 



21og(£>) 



l0 s((^T 



But now 



21og(£>) 



41og(D) 



lo g((^fcl 



as long as -h < 1 



< 



21og(d)-log(4)-log(d-l) 
AK \og(D) 



log(d) 



h . lo |^ , which finishes the proof of the 
Theorem. □ 



Appendix D: Some results for non translational 
invariant MPS 



Lemma 14. Let Ai,A 6 A4d, then there exist Bi,Ci S 
A4 d such that if we consider the state 



— ^ ^ r (^u ■ • • Ai L B iL+1 Ci L+2 ...C. 



»i ,...,»£ 



i N ) \H,---lN) 



(Dl) 



then the normalized reduced density matrix for L particles 
(particles 1-L) is 



Pi... 



L= J2 tr(At L ...A] i AA il ...A il ^\i 1 ,...i N ){j 1 ,...j L \ 

(D2) 



»1 ,...,»£ 



Proof. We consider the channel defined as 



l{X) = Y,ViXVl 



(D3) 



where V\\D is a diagonal unitary matrix with differ- 
ent incommensurable eigenvalues (such that Vf still has 



10 



different eigenvalues for all k € N), V2V-D is a random 
unitary matrix with non-zero entries and Vi — 0d,i € 
{3, This channel is trace preserving and unital. 

On the one hand, it is trivial to see that the only matri- 
ces that commute with V\ are diagonal matrices. On the 
other hand, to find which of these diagonal matrices com- 
mute with V2 it is enough to consider the algebraic system 
of equations in coordinates for [V^-X] = from where 
we get that, since (Vz)^ ^ 0, and (X)ij = if i 7^ j, 
then (X)u — (X)jj = for all i ^ j. From this, we de- 
duce that the only matrices that commute with all of 
the Kraus operators for our channel are multiples of the 
identity matrix. Liidcrs' Theorem [25] guarantees that 
our channel has the identity as its unique fixed point. 
Since E is an irreducible channel [26], all its eigenvalues 
of modulus 1 are /c-roots of unity, where k £ {1, . . . ,D 2 }. 
Let Y be such that E(Y) = aY for \a\ = 1. It is 
clear that E k (Y) = Y and, again by Liiders' Theorem, 



[Vf,Y] = = [V 2 Vf-\Y]. Reasoning as above, Y is a 
multiple of the identity, which implies that a — 1; hence, 
the channel is primitive [26]. 
We can define now 



B 4 



Cfc = Vk 



(D4) 



where Vi are the Kraus operators for our channel. If we 
consider the state 

\i>)= E te(A il ...A iL B iL+1 C iL+2 ...C iN )\i 1 ,...i N ) 



tl,...,»L 

ii+i, -W 



(D5) 

and compute the normalized reduced density matrix for 
particles 1...L, we obtain 



PL 



E 



J2(a\ [a h ...A 1l E b E, 



n-(l+v 
c 



(\a){p\)A] L ---A] 1 



11, 



■ ,IL 

•Jl 



It is clear that E^ ^ L+1 \\a) (f3\) = 5 a pl, up to an expo- nentially small correction. This leads us to 



Pl= £ l^\[A il ...A iL E B (l)A] L ...A] i ]\f3)\ 

ii,...,i L \a,/3 ) 

= E Enh,.-4A4...4]w 

ii,...,iL \a,/3 J 

once again, up to a exponentially small correction, just 
as we wanted to prove. □ 
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